Buscando planetas

En este notebook jugaremos con los niveles de una imagen reducida de Fomalhaut para encontrar al exoplaneta Fomalhaut b. El archivo .fits de imagen utilizado en este notebook es PCA_FOMAL_full_centroids.fits esta bajo embargo, pero si se necesita para propositos educativos, por favor envie un correo a lhquirogan@gmail.com solicitando este archivo.



In [1]:

    
from astropy.io import fits
import matplotlib.pyplot as plt
%matplotlib inline



In [2]:

    
n=fits.open('PCA_FOMAL_full_centroids.fits')



In [3]:

    
dst=n[0].data



In [4]:

    
dst









    Out[4]:





array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
       [ 0.,  0.,  0., ...,  0.,  0.,  0.],
       [ 0.,  0.,  0., ...,  0.,  0.,  0.],
       ..., 
       [ 0.,  0.,  0., ...,  0.,  0.,  0.],
       [ 0.,  0.,  0., ...,  0.,  0.,  0.],
       [ 0.,  0.,  0., ...,  0.,  0.,  0.]], dtype=float32)



In [5]:

    
plt.imshow(dst)









    Out[5]:





<matplotlib.image.AxesImage at 0x7fa952e52f60>



In [16]:

    
dst.min()









    Out[16]:





-0.70390803



In [15]:

    
dst.max()









    Out[15]:





0.96722591



In [13]:

    
ninf=0.
nsup=0.005



In [14]:

    
plt.figure(figsize=(10,10))
plt.imshow(dst*(dst>ninf)*(dst<nsup),cmap="gist_heat")
plt.scatter(1173,1205,s=300,facecolors='none', edgecolors='b',lw=2)
plt.xlim(900,1250)
plt.ylim(1100,1400)









    Out[14]:





(1100, 1400)



In [6]:

    
plt.figure(figsize=(10,10))
plt.imshow(dst*(dst>-0.0045)*(dst<0.004),cmap="gist_heat")
plt.scatter(1173,1205,s=300,facecolors='none', edgecolors='b',lw=2)
plt.xlim(900,1250)
plt.ylim(1100,1400)









    Out[6]:





(1100, 1400)



In [ ]: